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1. INTRODUCTION 

This study is related to the fluctuation theory 
of electromagnetic (EM) fields, charges, and currents 
in systems formulated in the three-dimensional (3D) 
'.Cartesian space of points r — (x,y,z). The studied 
problem is inhomogeneous along one of the axis, say 
along the first coordinate x, and translationally invari- 
ant along the plane formed by the remaining two co- 
ordinates normal to x, denoted as R = {y,z). The 
model consists of two semi-infinite media (conductor, 
■ dielectric or vacuum) with the frequency-dependent di- 
'electric functions ei(w) and e2{'^) which are localized 
,in the complementary half-spaces a; > and x < 0, 
'respectively (see Fig. 1). The interface between media 
is the plane x — Q. We assume that the media have 
no magnetic structure and the magnetic permeabili- 
,ties fJ,i = IJ,2 — 1- The media and the radiated EM 
'field are in thermal equilibrium at some temperature 
T, or the inverse temperature /? = 1/(^3^) with fc^ 
being the Boltzmann constant. 
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1. Two semi-infinite media characterized by the 
dielectric functions ei{uj) and e2(w). 
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The different electric properties of the media give 
.rise to a surface charge density a{t, R) at time t and at 
a point r — (0,R) on the interface which must be un- 
.derstood as being the microscopic volume charge den- 
'sity integrated along the a;- axis on some microscopic 
depth from the interface. In thermal equilibrium, this 
'quantity fluctuates in time around its mean value, usu- 
ally equal to 0. Based on the elementary electrodynam- 
ics, the surface charge density is associated with the 
discontinuity of the normal x component of the elec- 
tric field at the interface: In Gauss units, we have 



47rcr(t,R) = £;+(t,R) --B^(t,R), 



(1) 



interface through the limit x 0+ {x 0~). The 
tangential y and z components of the electric field are 
continuous at the interface. 

The fluctuations of the surface charge density at 
two points on the interface, with times different by t 
and distances different by i? = |R|, are correlated due 
to the EM interactions among the charged particles 
forming the two media. We are interested in the (sym- 
metrized) two-point surface charge correlation function 
defined by 



5(t,R)= i(a(t,R)a(0,0) 



a{0,0)a{t,R)f, (2) 



where (• • •)"^ represents a truncated equilibrium sta- 
tistical average, (AB)'^ = (AB) - {A){B), at the in- 
verse temperature (3. With regard to the relation ([T]), 
the function ([2]) is related to fluctuations of the EM 
field on the interface between media. Although the sys- 
tem is not in the critical state, due to the presence of 
long-ranged EM forces in the media the asymptotic 
large-distance behavior of the surface charge correla- 
tion function ([2]) exhibits in general a long-ranged tail 
of type 

Mt) 

i?3 ' 



/3S{t, R) 



R 



(3) 



where the form of the prefactor function h{t) depends 
on "the physical model" used. We can consider: 

• Either classical or quantum mechanics. Accord- 
ing to the correspondence principle, the two 
theories provide the same results in the high- 
temperature limit f3h — !■ 0. 

• Either non- retarded or retarded regime of the EM 
interaction. In the non-retarded (non-relativistic) 
regime, the speed of light c is taken to be in- 
finitely large, c — oo, ignoring in this way mag- 
netic forces acting on the particles. In the re- 
tarded (relativistic) regime, c is assumed finite 
and so the particles are fully coupled to both 
electric and magnetic parts of the radiated EM 
field. 

The particles themselves are non-relativistic, classical 
or quantum. 

It is useful to introduce the Fourier transform 



S{t,q)= Jd^R exp(-iq-R)S'(i,R), 



(4) 



where the superscript -I- (— ) means approaching the 



with q = {qy,qz) being a 2D wave vector. Since, in 
the sense of distributions, the 2D Fourier transform of 
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\/R^ is — 27rg, a result equivalent to ([3]) is that /3S'(t, q) 
has a kink singularity at q = 0, behaving like 



l3S{t,q) 2TTh{t)q, 



0. 



(5) 



The correlation function S{t, R) and its long- 
ranged decay ^ are of importance for two reasons. 
Firstly, the Fourier transform is a measurable quantity 
by scattering experiments. The fact that S{t, q) is lin- 
ear in q for small q makes this quantity very different 
from the usual correlation functions with short-ranged 
decay which are proportional to q^ in the limit q — > 0. 
The second application arises in the context of the av- 
erage polarization P of a charged system, localized in- 
side a domain A of volume | A| entoured say by vacuum 
and exposed to a constant electric field E. In the linear 
limit, we have 



|A| 



E 

k 



XjkEk 



J, 



x,y, 



(6) 



where Xjk is the dielectric susceptibility tensor. One 
expects that an extensive quantity like P divided by 
the volume |A| tends, for a fixed temperature, to a 
constant in the thermodynamic limit |A| oo, inde- 
pendently of the shape of the infinite domain A. This is 
not true in the case of the dielectric susceptibility ten- 
sor. It can be shown that Xjk contains contributions 
from both short-ranged bulk and long-ranged surface 
charge ^ correlation functions which causes its shape 
dependence [TJ H] . Such a phenomenon is predicted by 
the macroscopic laws of electrostatics [3]. 

The two media configuration studied so far was re- 
stricted to a conductor, localized say in the half-space 
X > with the dielectric function ei{uj) = e{uj), in con- 
tact with vacuum of the dielectric constant €2(^0) ~ 1. 
The obtained results can be summarized as follows. 

• Classical non-retarded regime: Let the con- 
ductor be modelled by a classical Coulomb fiuid com- 
posed of charged particles with only the instantaneous 
Coulomb interactions. By a microscopic analysis [4J, 
the long-ranged decay of the static (i.e. t = 0) surface 
charge correlation was found such that 



1 

'8^' 



(7) 



where the subscript "cl" means "classical" and the up- 
perscript "nr" means "non-retarded" . The same result 
has been obtained later by simple macroscopic ar- 
guments based on a combination of the linear response 
theory and the electrostatic method of images. Note 
the universal form of ft.^"'^''(0), independent of the com- 
position of the classical Coulomb fluid. 

• Quantum non-retarded regime: The quan- 
tum description of a general conductor is very com- 
plicated. A simplification arises in the so-called jel- 
lium model, i.e. a system of identical pointlike par- 
ticles of charge e, mass m and bulk number density 



n, immersed in a uniform neutralizing background of 
charge density —en. The dynamical properties of the 
jellium have a special feature: There is no viscous 
damping of the long-wavelength plasma oscillations 
for identically charged particles [6] . The frequencies of 
non-retarded nondispersive long-wavelength collective 
modes, namely lo^ of the bulk plasmons and lOg of the 
surface plasmons, are given by 



/ 47rne^ 
\ m 



1/2 



(8) 



The dielectric function of the jellium is well described 
by a simple one-resonance Drude formula 



e(c^) = 1 - 



ijj{uj -f i?7) ' 



(9) 



where the dissipation constant 77 is taken as positive 
infinitesimal, 77 — > 0"*". The prefactor function /i(t) for 
the quantum jellium was obtained in Refs. [71IB]. It has 
the nonuniversal form 



where 



[2(7(ti;s) cos(wsi) — g(y-p) cos(wpt)] 



(10) 
(11) 



and the subscript "qu" means "quantum" . In the high- 
temperature limit (5% — > 0, the function g{u}) = 1 for 
any uj and the quantum formula (jlOp reduces to the 
classical non-retarded one 



^ir'(0 = [2cos(wsi) - cos(tjj,t)] 



1 



(12) 



For t = 0, we recover the classical static result ©. 

The present study reviews some new results in the 
field Uni [H]- These results concern a generahza- 
tion of the formalism to contacts between all kinds of 
media, i.e. conductors, dielectrics and vacuum, and to 
the retarded regime. It turns out that, for any time 
the inclusion of retardation effects surprisingly con- 
verts the quantum results into the static classical ones. 

To deal with the physical problem of such complex- 
ity, we use a macroscopic theory of equilibrium ther- 
mal fluctuations of EM field, initiated by Rytov [12] 
and further developed in Ref. [13]; we adopt the nota- 
tion from the course by Landau and Lifshitz [T3] (see 
Section 2). Based on Rytov's fluctuational theory, we 
derive in Section 3 the general formula for the surface 
charge correlation function between two semi-infinite 
media. The analysis of the asymptotic large-distance 
form of this formula, in both non-retarded and retarded 
regimes, is the subject of Section 4. We first consider 
the special configuration of the jellium conductor in 
contact with vacuum (Subsection 4.1) and then the 
general configuration (Subsection 4.2). Section 5 is a 
Conclusion. 
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2. FLUCTUATIONAL ELECTRODYNAMICS 

Let us consider a non-magnetoactive {fJ- — 1) 
medium with the isotropic, frequency and (possibly) 
position dependent, dielectric function e{uj;r). The 
coupled EM field is defined, in the classical format, 
by the scalar potential (j>{t, r) and the vector poten- 
tial A{t,r). In the considered Weyl gauge = 0, the 
microscopic electric and magnetic fields are given by 



E 



1 dA 

c~dt' 



B = curl A. 



(13) 



The elementary excitations of the quantized EM field 
are described by the photon operators Aj (j = x, y, z) 
which are self-conjugate Bose operators; Aj{t,Y) will 
denote the vector-potential operator in the Heisenberg 
picture. The construction of all types of photon Green's 
functions is based on the retarded Green function D, 
defined as the tensor 

^i?,fc(^;r,r') = (i,(t, r)ifc(0, r') - ife(0, r')i,(i, r)) 

(14) 

for t > and equal to for t < 0. The Fourier trans- 
form in time of the retarded Green function reads 



/"OO 

i^,fe(w;r,r') = / e''"*-D,fe(t; r, r'). 
Jo 



(15) 



For media with no magnetic structure, the Green func- 
tion tensor possesses the important symmetry 



Djk{(^;Y,Y') = Dkj{uj;r' ,r) 



(16) 



The EM fields, being in thermal equilibrium with 
medium, are random variables which fluctuate around 
their mean values obeying macroscopic Maxwell's 
equations. To deal with thermal fluctuations, Rytov 
put a classical current j(i,r) due to the thermal mo- 
tion of particles into the medium [12j . This current acts 
as an "external force" on the vector-potential operator 
in the interaction Hamiltonian 

Hint(t) = -i j d\i{t,r)-k. (17) 

The mean values of the components of the vector- 
potential operator can be expressed in terms of the re- 
tarded Green function by using Kubo's linear response 
in currents: 



^ I. 



(18) 

The fact that the mean value A satisfies the macro- 
scopic Maxwell equation due to the classical current 
j implies the differential equation of the dyadic type 
fulfilled by the retarded Green function tensor: 



3 

E 

1=1 



dxj dxi 



SjiA- 6ji—e{uj;r) 



Dik{uj]Y,v') 



Here, in order to simplify the notation, the vector 
r — [x, y, z) is represented as (xi,X2,x^). The differen- 
tial equation must be supplemented by certain bound- 
ary conditions. The second space variable r' and the 
second index k only act as parameters, the boundary 
conditions are formulated with respect to the coordi- 
nate r and the Green function Dik{uj;r,r') is consid- 
ered as a vector with the components I = x,y, z. There 
is an obvious boundary condition of regularity at in- 
finity, |r| — > oo. At an interface between two differ- 
ent media, the boundary conditions correspond to the 
macroscopic requirements that the tangential compo- 
nents of the fields E and H = B be continuous. The 
vector potential is related to the Green function tensor 
via the linear response ([18]), and the electric and mag- 
netic fields are related to the vector potential by 
This is why the role of the vector components Ei , up 
to an irrelevant multiplicative constant, is played by 
the quantity 



-D(fc(w;r,r') 



(20) 



and the role of the vector component Hi is played by 
the quantity 



]cmlijDjk{^^;r, r'). 



(21) 



Here, we use the notation curlj^ = Y^^eimjd/dxm 
with eimj being the unit antisymmetric pseudo-tensor. 
For our geometry in Fig. 1, the tangential components 
(|20|) and (|2T|) . which are continuous at the interface 
X — correspond to indices I = y, z. 

Eq. (fT5)) defines —Djk{uj;r,r')/{hc'^) as the ten- 
sor of generalized susceptibilities corresponding to the 
variable A/c. The fluctuation-dissipation theorem tells 
us that the fluctuations of random variables can be ex- 
pressed in terms of the corresponding susceptibilities. 
For the assumed symmetry (I16D . the theorem implies 
that 

[Aj{r)Ak{r')]^ = - coth(/3?ia;/2) Iml?,fc(cj; r, r'), 

(22) 

where the spectral distribution [Aj{r)Ak{r')]uj is the 
Fourier transform in time of the symmetrized (trun- 
cated) correlation function 

i (i,(<,r)ifc(0,r') +ifc(0,r')i,(t,r))^ . (23) 

The spectral distribution of the electric-field fluctua- 
tions can be easily found by using the first relation in 
Eq. ([I3D, 



[A,{r)A,{r% 



r.2 L^-^J 



-47rMjfc(5(r- r'). (19) 



-coth(/3?icj/2) 
xlmi^,fe(w;r,r'). (24) 
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3. SURFACE CHARGE CORRELATIONS 

The retarded Green function tensor for the studied 
problem of two semi- infinite media in Fig. 1 is obtained 
as the solution of the dyadic difi^erential Eq. ([19]) sup- 
plemented by the mentioned boundary conditions [9^. 
Since the system is translationally invariant in the R- 
plane perpendicular to the x axis, we introduce the 
Fourier transform of the retarded Green function ten- 
sor with the wave vector q — {qy, q^), 



„iq(R-R' 



Djk{uj,c\:,x,x'). 



(25) 

Let us define for each of the half-space regions the in- 
verse length Kj{uj,q) (j = 1,2) by 



K'^j{uj,q) =q'^ -^ejiuj), Re Kj{lj, q) > 0; (26) 



from two possible solutions for Kj we choose the one 
with the positive real part in order to ensure the reg- 
ularity of the Green function at asymptotically large 
distances from the interface. We shall only need the 
quantity D^xii-^, Q, x, x') for which the obtained results 
can be summarized as follows: 
(i) lix.x' > 0, 



do not play any role in the large- distance asymptotic. 
Since the combination 



4:nh{cqf 1 



£2 _|_ £1 

we finally arrive at the quantum result 



- 2 



/3S'qu(t, q) 



e-'-*Im/H, 



(30) 



(31) 



where, in the retarded regime, f{uj) = /qu with 



/■II- J 

J qu 



(r) 



47r2 Ki{uj,q)e2{uj) + K2{LU,q)ei{uj) 



ei(w)e2(w) 



(32) 



In the non-retarded case, f{to) = /qS''' is obtained from 
(|32p by setting the speed of light c 00. According to 
Eq. (|26p . the inverse lengths ki = K2 = q in this limit, 
so that 



1 



1 



ei(w) e2(t^) ei(tj) + e2(w) 



(33) 



■6{x-x') 



2iTfi{cqY 



4. ANALYSIS OF THE RESULTS 



Kie2J_K2ei^_Ki(2;+a;') 
Kie2 + K2ei 



(ii) If X < and x' > 0, 

47r?i(cg)2 1 



^ Kie2 -I- K2ei 



(27) 



(28) 



(iii) The case x > and a;' < is deducible from 
Eq. (pS)) by using the symmetry relation 

(iv) 11 x^x' < 0, considering the 1^2 media ex- 
change symmetry, we obtain from Eq. (|27p that 



w^e2 



■5{x-x') 



UJ^e2K2 



K2€l - Kie2 ^^^(x+x') 
K2t\ -I- Kie2 



(29) 



The symmetrized surface charge correlation func- 
tion ([2]) is expressible in terms of the fluctuations of 
the symmetrized electric field a;a;-components by using 
relation ^ . These electric- field fiuctuations are related 
to the XX elements of the Green function tensor via Eq. 

The terms proportional to 5(x - x') in Eqs. (P7)) 
and (|29p can be ignored since they originate from the 
short-distance terms proportional to 5{r — r') which 



4-.1- Jellium conductor in vacuum 

We first considered the previously studied configu- 
ration of the jellium conductor of the dielectric func- 
tion ei((jj) given by the Drude formula ^ in contact 
with the vacuum of e2(w) = 1 [TU]. 

In the non-retarded regime, using the Weierstrass 
theorem 



lim 

r/^0+ X ± 177 



= r 



^ mS{x) 



(34) 



{V denotes the Cauchy principal value) in the represen- 
tation the integration over frequency (Plj) leads 
the previous result pH)) . This result is valid for inter- 
mediate distances R on the interface, given by the in- 
equalities Aph <C R ^ c/o^p, where Aph oc /She stands 
for the thermal de Broglie wavelength of photon and 
c/u>p is the wavelength of electromagnetic waves emit- 
ted by charge oscillations at frequency tOp. 

The analysis of the frequency integral ([51]) with the 
retarded function (|32[) is much more complicated [3]. 
At specific values of to, given by the surface-plasmon 
dispersion relation 



g^(ei + £2) 5-eie2 = 0, 



(35) 
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the integration passes in the infinitesimal vicinity of a 
singularity. The final result for the small (/-expansion 
is of the classical static type 



/3S«(t,(7) = 



47r 



o(q). 



(36) 



In other words, for both static and time-dependent sur- 
face charge correlation functions, the inclusion of retar- 
dation effects causes the quantum prefactor to take its 
universal static classical form 



1 



(37) 



This formula is valid in the strictly asymptotic region 
of distances R on the interface, given by the inequality 
c/ujp <C R- Note that it holds 



/^i;^(o) = /^ir^(o)^/ici(o). 



(38) 



This relation is of general validity due to the Bohr-van 
Leeuwen theorem [151 H] • 

The same result ([57)) was derived by another (par- 
tially microscopic) method, based on the analysis of 
the collective vibration modes of the system [TO] . 

4-2. General configuration 

For an arbitrary configuration of the plane contact 
between two distinct media, the small-q analysis of the 
surface charge correlation function (pijl was done in 
Ref. [TT] by using general analytic properties of di- 
electric functions in the complex frequency upper half- 
plane and contour integration techniques. 

We start with the simpler static case t = 0. The 
formula ([5T|) reads 



/35qu(0,g) 



dw 



Im/(a;) 



(39) 



where the function /(w) is given by (l32l) in the retarded 
regime and by p3p in the non-retarded regime. For 
a real frequency uj, the symmetry relations e*(w) — 
e(— w) and k*{uj) — k{—uj) imply f*{oj) = f{~uj), i.e. 



Re/(L^) = Re/(-a;), Im /(w) - -Im /(-c^). 



As w ^ 0, Im/(0) = and 



(40) 



Re/(0) = 



47r2 



1 



1 



ei(0) £2(0) ei(0) + e2(cc;) 



(41) 

Let Ci be the path following the real axis, except 
it goes around the origin w = in a small semicircle in 
complex upper half-plane whose radius e tends to zero 
and C2 be the closure of Ci by a semicircle at infinity 
(see Fig. 2). The integration along the real axis in Eq. 



Im CO 




Re CO 



Fig. 2. The contour in the complex frequency plane for 
t = 0. 



([5^ is expressible in terms of the path integration over 
the closed contour C — Ci U C2 as follows 



Im/(c^) = 71/(0) 



Im (f — /(w) 
Jc ^ 



(42) 



The integral over the contour C can be evaluated by 
using the residue theorem at poles {iOj} of the function 
f{uj) in the uj upper half-plane bounded by C, 



Im 



c 



dcj „, , „ v-^ 



(43) 



provided that 'R,es{f , uj j) / u) j is real (which will be the 
case); Res denotes the residue. Both the retarded (j32p 
and non-retarded (|33|) versions of the /-function con- 
tain g{uj) defined in Ijlip . Since g{uj) can be expanded 
in Lo as |17) 



5(c.) = l + ^- 



2tt ^ 



(44) 



it has in the upper half-plane an infinite sequence of 
simple poles at the imaginary Matsubara frequencies 



i^j, Res(g,Wj) 



(j = l,2,...). (45) 



Using general analytic properties of dielectric functions 
it was shown in Ref. [llj that, in both retarded and 
non-retarded regimes, these are the only poles of the 
/-function inside the contour C. The static correlation 
function ([55)1 is therefore expressible as follows 



/3^qu(0,g) 



q_ 

An 



£2(0) 



F(0,g) = 27r 



^i(O) 
Res(/, i^j 



£l(0) + £2(0) 

(46) 
(47) 



The first term on the rhs of Eq. ([46[l . linear in q, is 
independent of Ph and c, the g-dependence of the static 
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function F{0, q) depends on the considered (retarded 
or non-retarded) regime. 

In the retarded case ([5^ . we have 



27^^ '*i(*Cj>2(iCi) 
ei(jCj>2(iG-) 



K2(i?j)ei(i0) 
(48) 



We are interested in the Hmit q 

1/2 



for which 

Ki,2ii^]) ~ 0<^i!2'(«Cj) (ei,2(«^i) are real). Since oc j 
and — 1 = 0(l/j^) in the hmit j — > oo, the 

sum in (|48p converges. This means that the function 
Fqu (0, g), being of the order 0{q^), becomes neghgi- 
ble in comparison with the first term in Eq. (|46p when 
g ^ 0. We find the static /i-prefactor associated with 
the asymptotic decay to be 



1 



1 



ei(0) 62(0) ei(0) + e2(0) 



(49) 

Since this expression does not depend on the temper- 
ature and its classical PTi ^ limit is the same, i.e. 



,(nr) 



(0) = /lcl(O) with 



^cl(0)=-^ 



87r2 [ei(0) £2(0) ei(0) + e2(0) 

;(50) 

This classical result was derived independently in the 
non-retarded regime by using the linear response the- 
ory combined with the electrostatic method of images 
We see that from the dielectric functions of the 
two media only their values at zero frequency appear; 
we recall that e(0) ioo for conductors, e(0) = eo > 1 
for dielectrics and e(0) = 1 for vacuum. The /i-prefactor 
is nonzero for an arbitrary configuration of different 
media, except for the special case of two conductors. 
In the non- retarded case ([55)1 . we have 



oo r 



1 



1 
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e2(i^j) 



(51) 



ei(i^j) + e2(j0) 

It is evident from the asymptotic behavior e(z^j) — \ — 
0{l/j^) that the sum in ([5T|l converges. The function 



-Fqu'^''(0, (?) is thus of the order 0{q) and h'i^^O) is a 
complicated function of the temperature. 

A more laborious analysis is needed for the time- 
dependent case t llj. The final result is 



(nr) 



1 



1 



ei(0) 62(0) ei(0) + e2(0) 



(52) 

As before, the quantum time-dependent prefactor to 
the asymptotic decay takes, for any temperature, its 
classical static form ([50]) . 



5. CONCLUSION 

Although the present study is rather technical, its 
main result (|52p is of importance. We do not under- 
stand which are physical reasons that the inclusion of 
retardation effects causes the asymptotic decay of time- 
dependent quantum surface charge correlation func- 
tions to take its static classical form, independent of 
the temperature factor ph and the speed of light c. 

This is one of the rare phenomena in the Con- 
densed Matter physics when retardation (relativistic) 
effects play an essential role. 
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